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Abstract. Absorption imaging with quasi-resonant laser light is a commonly used 
technique to probe ultra-cold atomic gases in various geometries. Here we investigate 
some non-trivial aspects of this method when it is applied to in situ diagnosis of a 
quasi two-dimensional gas. Using Monte Carlo simulations wc study the modification 
of the absorption cross-section of a photon when it undergoes multiple scattering in 
the gas. We determine the variations of the optical density with various parameters, 
such as the detuning of the light from the atomic resonance and the thickness of the 
gas. Wc compare our results to the known three-dimensional result (Beer-Lambert 
law) and outline the specific features of the two-dimensional case. 
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1. Introduction 

The study of cold atomic gases has recently shed new light on several aspects of quantum 
many-body physics [H El El H] . Most of the measurements in this field of research are 
based on the determination of the spatial density of the gas [3] . For instance one can use 
the in situ steady-state atomic distribution in a trapping potential to infer the equation 
of state of the homogenous gas [6]. Another example is the time-of-fiight method, in 
which one measures the spatial density after switching off the trapping potential and 
allowing for a certain time of ballistic expansion. This gives access to the momentum 
distribution of the gas, and to the conversion of interaction energy into kinetic energy 
at the moment of the potential switch-off. 

To access the atomic density n{r), one usually relies on the interaction of the atoms 
with quasi-resonant laser light. The most common method is absorption imaging, in 
which the shadow imprinted by the cloud on a low intensity probe beam is imaged 
on a camera. The simplest modelling of absorption imaging is based on a mean-field 
approach, in which one assumes that the local value of the electric field driving an 
atomic dipole at a given location depends only on the average density of scatterers. 
One can then relate the attenuation of the laser beam to the column atomic density 
y) = J n{r) dz along the line-of-sight z. The optical density of the cloud 
V{x,y) = ln[/in(x, ?/)//out(a^, y)] is given by the Beer-Lambert law 

VBLix,y) = an^''"'\x,y), (1) 

where a is the photon scattering cross-section, and /in (resp. /out) are the incoming (resp. 
outgoing) intensity of the probe laser in the plane xy perpendicular to the propagation 
axis. For a closed two-level atomic transition of frequency Uq = ckQ, o depends on the 
wavelength Aq = 2tt /kQ associated to this transition and on the detuning A = w — 
between the probe light frequency u and the atomic frequency: 

Here F represents the natural line width of the transition {i.e., is the natural life 
time of the excited state of the transition). Eq. ([2]) assumes that the intensity of the 
probe beam is much lower than the saturation intensity of the atomic transition. Quasi- 
resonant absorption imaging is widely used to measure the spatial distribution of atomic 
gases after a long time-of-fiight, when the density has dropped sufficiently so that the 
mean-field approximation leading to Eq. ([T]) is valid. 

One can also use absorption imaging to probe in situ samples, at least in the case 
where an^^"^^ is not very large so that the output intensity is not vanishingly small. 
This is in particular the case for low dimensional gases. Consider for example a 2D 
gas, such that the translational degree of freedom along z has been frozen. For a probe 
beam propagating along this axis, one can transpose the Beer-Lambert law of Eq. ([T]) 
by simply replacing the column density by the surface density n^"^^^ of the gas. This 2D 
Beer-Lambert law can be heuristically justified by treating each atom as a disk of area a 
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that blocks every photon incident on it. In an area a containing N = An^'^^^ ^ 1 
randomly placed atoms, the probability that a photon is not blocked by any of the disks 
is (l-aM)^^exp(-an(2D)). 

In a quasi-2D gas there is however an important limitation on the optical densities to 
which one may apply the Beer-Lambert prediction of Eq. ([l]). Already for aQU^'^^^ = 1 
the mean interparticle distance is only 0.7 Aq and one may expect that the optical 
response of an atom strongly depends on the precise location of its neighbours. More 
precisely the exchange of photons between closely spaced atoms induces a resonant 
van der Waals interaction that significantly shifts the atomic resonance frequency with 
respect to its bare value Uq. The optical density of the gas at resonance may then be 
reduced with respect to Eq. Q, and this was indeed observed in a series of experiments 
performed with a degenerate ^^Rb gas [71 |H] ■ 

The general subject of the propagation of a light wave in a dense atomic sample, 
where multiple scattering plays an essential role, has been the subject of numerous 
experimental and theoretical works (see e.g. [HI in the context of cold atoms, and 
[TT] for a review). Here we present a quantitative treatment of the collective effects that 
appear when a weak probe beam interacts with a quasi-2D atomic gas. We consider an 
ensemble of N atoms at rest with random positions and we investigate the transmission 
of quasi-resonant light by the atom sheet. We model the resonance transition between 
the atomic ground (g) and excited (e) states by a J^, = o Jg = 1 transition. We present 
two equivalent approaches; the first one is based on the calculation of the field radiated 
by an assembly of N dipoles, where each dipole is driven by an external field plus the field 
radiated by the A^— 1 other dipoles; the second one uses the standard T matrix formalism 
of scattering theory. We show that in both cases the optical density of the medium can 
be determined by solving the same 3A^ x 3A^ linear system. A similar formalism has 
been previously used for the study of light propagation in small 3D atomic samples, in 
the presence of multiple scattering (see e.g. [12], [131 El CSl UHl HTJ HH HE]). However 
its application to quasi-2D samples has (to our knowledge) not yet been investigated, 
except in the context of Anderson localisation of light [13] . Our numerical calculations 
are performed for A^ = 2048 atoms, which is sufficient to reach the 'thermodynamic limit' 
for the range of parameters that is relevant for experiments. We show in particular that 
even for moderate values of aon^^°^\ the optical density is notably reduced compared 
to what is expected from the Beer-Lambert law {e.g., more than 20% reduction for 
cTon'-™'^ = 1). We investigate how the absorption line shape is modified by the resonant 
van der Waals interactions and we also show how the result ([T]) is recovered when one 
increases the thickness of the gas, for a given column density n^'^"^\ 

The paper is organised as follows. In section |2j we detail the modelling of the atom- 
light interaction with the two-level and rotating wave approximations. Then we explain 
the principle of the calculation for the absorption of a weak probe beam crossing the 
atom slab (section |3]). The ensemble of our numerical results are presented in section |4j 
Finally in section |5] we discuss some limitations to our model and draw some concluding 
remarks. 
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2. Modelling the atom-light interaction 

2.1. The electromagnetic field 

We use the standard description of the quantised electromagnetic field in the Coulomb 
gauge [20], and choose periodic boundary conditions in the cubic-shaped quantisation 
volume V = L^LyLz. We denote aq s the destruction operator of a photon with wave 
vector q and polarisation s (s±qr). The Hamiltonian of the quantised field is 

Hf = ^ hcqa\^^aq^s , (3) 

q,s 

and the transverse electric field operator reads E{r) = E^^\r) + E^^\r) with 

q,s 

and E^-\r) = (^E'-+\r)^ . The wave vectors q are quantised in the volume V as 
Qi = 2Trni/Li, i = x,y, z, where Ui is a positive or negative integer. 

2.2. The atomic medium 

We consider a collection of N identical atoms at rest in positions Vj, j = 1, . . . N. We 
model the atomic resonance transition by a two- level system with a ground state \g) 
with angular momentum Jg = and an excited level of angular momentum Je = 1. We 
choose as a basis set for the excited manifold the three Zeeman sublevels le^), a = x,y, z, 
where \ea) is the eigenstate with eigenvalue of the component Ja of the atomic angular 
momentum operator. We denote hwQ the energy difference between e and g. The atomic 
Hamiltonian is thus (up to a constant) 

N 

= 5Z ^o\j ■ ea){j ■■ e^l- (5) 

j=l a=x,y,z 

The restriction to a two-level approximation is legitimate if the detuning A between 
the probe and the atomic frequencies is much smaller than uq. The modelling of this 
transition by a Jg = Je = 1 transition leads to a relatively simple algebra. The 
transitions that are used for absorption imaging in real experiments often involve more 
Zeeman states ( J^, = 2 ^ Jg = 3 for Rb atoms in [HIB]), but are more complex to handle 
[21] 122] and they are thus out of the scope of this paper. However we believe that the 
most salient features of multiple scattering and resonant Van der Waals interactions are 
captured by our simple level scheme. 

2.3. The atom-light coupling 

We treat the atom-light interaction using the electric dipole approximation (length 
gauge), which is legitimate since the resonance wavelength of the atoms Aq is much 
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larger than the atomic size. We write the atom-hght couphng as: 

V = -Y,DrE{r,), (6) 

j 

where Dj is the dipole operator for the atom j. We will use the rotating wave 
approximation (RWA), which consists in keeping only the resonant terms in the coupling: 

V^-Y, D^^^ ■ E^^^ i^j) + h-c. , (7) 
j 

where h.c. stands for Hermitian conjugate. Here D^^^ represents the raising part of the 
dipole operator for atom j: 

^r^=^ E \j--ea){j:9\uc., (8) 

a=x,y,z 

where d is the electric dipole associated to the g — e transition and iia is a unit vector 
in the direction a. 

When a single atom is coupled to the electromagnetic field, this coupling results 
in the modification of the resonance frequency (Lamb shift) and in the fact that the 
excited state e acquires a non-zero width T 

(9) 



For simplicity we will incorporate the Lamb shift in the definition of Uq. Note that the 
proper calculation for this shift requires that one goes beyond the two-level and the 
rotating wave approximations. The linewidth F on the other hand can be calculated 
from the above expressions for V using the Fermi golden rule. 

The RWA provides a very significant simplification of the treatment of the atom- 
light coupling, in the sense that the total number of excitations is a conserved quantity. 
The annihilation (resp. creation) of a photon is always associated with the transition 
of one of the atoms from g to e (resp. from e to g). This would not be the case if 
the non-resonant terms of the electric dipole coupling d\'^'^ ■ E^~^ and d[ ^ ■ E''^^ were 
also taken into account. The small parameter associated to the RWA is A/uq, which is 
in practice in the range 10"^ — 10^^; the RWA is thus an excellent approximation. 

Formally the use of the electric dipole interaction implies to add to the Hamiltonian 
an additional contact term between the dipoles (see e.g. [231 112] )• This term will play 
no role in our numerical simulations because we will surround the position of each atom 
by a small excluded volume, which mimics the short range repulsive interaction between 
atoms. We checked that the results of our numerical calculations (see Sec. ^ do not 
depend on the size of the excluded volume, and we can safely omit the additional contact 
term in the present work. 

3. Interaction of a probe laser beam with a dense quasi-2D atomic sample 



We present in this section the general formalism that allows one to calculate the 
absorption of a quasi-resonant laser beam by a slab of atoms. We address this 
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question using two different approaches. The first one maps the problem onto the 
collective behaviour of an assembly of oscillating dipoles [12]. The equation of motion 
for each dipole is obtained using the Heisenberg picture for the Hamiltonian presented 
in section [2j It contains two driving terms, one from the incident probe field and one 
from the field radiated by all the other dipoles at the location of the dipole under study. 
The steady-state of this assembly of dipoles is obtained by solving a set of 3A^ linear 
equations. The second approach uses the standard quantum scattering theory [21], 
which is well suited for perturbative calculations and partial resummations of diagrams. 
We suppose that one photon is incident on the atomic medium and we use resummation 
techniques to take into account the multiple scattering events that can occur before 
the photon emerges from the medium. The relevant quantity in this approach is the 
probability amplitude Ta that the outgoing photon is detected in the same mode as the 
incident one [H] [IT] , and we show that Tu is obtained from the same set of equations 
as the values of the dipoles in the first approach. 



3.1. Wave propagation in an assembly of driven dipoles. 

In this section we assume that the incident field is prepared in a coherent state 
corresponding to a monochromatic plane wave Ei^ e e*('^^~'^*). We choose the polarization 
e to be linear and parallel to the x axis (e = Ux). Since we consider a Jg = ^ Je = 1 
transition, this choice does not play a significant role and we checked that we recover 
essentially the same results with a circular polarisation. Note that the situation would be 
different for an atomic transition with larger Jg and Jg since optical pumping processes 
would then depend crucially on the polarisation of the probe laser. 

The amplitude Ej^ is supposed to be small enough that the steady-state populations 
of the excited states Cj^a are small compared to unity. This ensures that the response of 
each atomic dipole is linear in this approximation is valid when the Rabi frequency 
dEi^/H is small compared to the natural width F or the detuning A. 

Using the atom-light coupling (|6|, the equations of motion for the annihilation 
operators aq s in the Heisenberg picture read: 



d,At) = cqa,,,{t) + Yl ■ ■ (10) 

' j' 

This equation can be integrated between the initial time to aiid the time t, and the 
result can be injected in the expression for the transverse field to provide its value at 
any point r: 

-t-to 



Eair,t) = Efrec,a(r,t) + 



cq 



j',a' q,s 



2eoV 

[tDj',a'it - t) e^'/-('--/)--3- s^sl, + h.c] , (11) 

where Efj-ee stands for the value obtained in the absence of atoms. We now take the 
quantum average of this set of equations. In the steady-state regime the expectation 
value of the dipole operator Dj(t) can be written dje~^'^^ + c.c, and the average of 
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Eireeif^t) IS the incident field e e*'-^^"'^*-* + c.c. . We denote the average value of the 
transverse field operator in r as {E{r, t)) = E{r) e"*"^* + c.c, and we obtain after some 
algebra (see e.g. [El [25]) 

Ea{r) = Ei^ e« e'''^ + 2_, 9a,a'{uj>) df^a' , (12) 

j',a' 

where we set uj = k{r — vj) (with k ^ k^ 



and 



0), 

9a,a'{u) = 6a,a'hl{u) + ^^^^/l2(u), (13) 

hi{u) = -^{u^ + iu-l), h2{u) = -%{-u^ -3iu + 3). (14) 

The function ga,a'{kr) is identical to the one appearing in classical electrodynamics [25], 
when calculating the field radiated in r by a dipole located at the origin. 

We proceed similarly for the equations of motion for the dipole operators D^j~^ and 
take their average value in steady-state. The result can be put in the form |i2j 



{S + i)dj^a+ ^ ga,a'{ujj')dj'^a' = -^^^ Ei.eae'''''' , (15) 

where the reduced detunig 5 = 2A/T has been defined in Eq. ^ and Ujj' = k{rji — Vj). 
This can be written with matrix notation 

[M]\X) = \Y) (16) 

where the 3A^ vectors \X) and \Y) are defined by 

Xj a 7^ ^ dj aj ^CfC ; (1''') 

evreo^L 

and where the complex symmetric matrix [M] has its diagonal coefficients equal to 6 + i 
and its off-diagonal coefficients (for j ^ j') given by ga,a'iujji). This matrix belongs to 
the general class of Euclidean matrices [27], for which the (i, j) element can be written 
as a function -F(rj, Vj) of points Vi in the Euclidean space. The spectral properties of 
these matrices for a random distribution of the r^'s (as it will be the case in this work, 
see Sec. |4|) have been studied in [271 EH [29l |30] . 



Eq. (15) has a simple physical interpretation: in steady-state each dipole dj is 
driven by the sum of the incident field E]^ and the field radiated by all the other dipoles. 
This set of 3A^ equations was first introduced by L. L. Foldy in [31] who named it. 



together with Eq. (12), "the fundamental equations of multiple scattering". Indeed for 



a given incident field, the solution of (16) provides the value of each dipole dj, which 



can then be injected in (12) to obtain the value of the total field at any point in space. 
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Figure 1. Two possible setups for measuring the absorption of an incident probe 
beam by a slab of atoms using a lens of focal /. a) Global probe, b) Local probe. 



3.2. Absorption signal 

From the expression of the average value of the dipoles we now extract the absorption 
coefficient of the probe beam and the optical density of the gas. We suppose that the 
atoms are uniformly spread in a cylinder of radius R along the z axis and located 
between z = —£/2 and z = i/2. We can consider two experimental setups to address this 
problem. The first one, represented in Fig. [T^, consists in measuring after the atomic 
sample the total light intensity with the same momentum k = ku^ as the incident probe 
beam. This can be achieved by placing a lens with the same size as the atomic sample, 
in the plane z = i' > i/2 just after the sample. The light field at the focal point of 
the lens F gives the desired attenuation coefficient. We refer to this method as 'global', 
since the field E{F) provides information over the whole atomic cloud. One can also use 
the setup sketched in Fig. which forms an image of the atom slab on a camera and 
provides a 'local' measurement of the absorption coefficient. In real experiments local 
measurements are often favored because trapped atomic sample are non homogeneous 
and it is desirable to access the spatial distribution of the particles. However for our 
geometry with a uniform density of scatterers, spatial information on the absorption 
of the probe beam is not relevant. Therefore we only present the formalism for global 
measurements, which is simpler to derive and leads to slightly more general expressions. 
We checked numerically that we obtained very similar results when we modelled the 
local procedure. 

We assume that the lens in Fig. [l^ operates in the paraxial regime, i.e., its focal 
length / is much larger than its radius R. We relate the field at the image focal point 
of the lens to the field in the plane z = i' just before the lens: 

^(^) = -TT / ^(^'?/'^') dxdy, (18) 

where the integral runs over the lens area. Since the incident probe beam is supposed to 
be linearly polarised along x, we calculate the x component of the field in F. Plugging 
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the value of the field given in Eqs. ( 12p!7 ) we obtain the transmission coefficient 



-^a; I, -r J I no atom j ^ -J 

This result can be simplified in the limit of a large lens by using an approximated 
value for the integral appearing in (19). We suppose that M' 1 so that the 
dominant part in g^^a is the e*"/M contribution to hi. More precisely the domain in 
the lens plane contributing to the integral for the dipole j is essentially a disk of radius 
a/A(£' — Zj) ~ y/Xl' centered on {xj,yj). When this small disk is entirely included in 
the lens aperture, i.e., the larger disk of radius R centered on x = ?/ = 0, we obtain 

l^g.AKr - rj)] dxdy ^ ^^5,,,,e^'=(^'-.) . (20) 



We use the result (20 ) for all atoms, which amounts to neglect edge effects for the dipoles 



located at the border of the lens, and we obtain: 

r = 1 - ^aon(™')n, (21) 

with n*^™'^ = N/ ttR^ and where the coefficient 11 is defined by 

n = ^E^^-.-^"^''^- (22) 

j 

This coefficient captures the whole physics of multiple scattering and resonant van der 
Waals interactions among the atoms. Indeed one takes into account all possible 
couplings between the dipoles when solving the 3A^ x 3A^ system [M]|X) = \Y). Once 
T is known the optical density is obtained from 

V = ln\r\'\ (23) 

As an example, consider the limit of a very sparse sample where multiple scattering 
does not play a significant role (cto'^''™'^ ^1). All non-diagonal matrix elements in [M] 
are then negligible and [M] is simply the identity matrix, times i + 6. Each Xj^^ solution 



of the system (16) is equal to /{i + 5), and we obtain as expected: 
3. 3. Light absorption as a quantum scattering process 

In order to study the attenuation of a weak probe beam propagating along the z axis 
when it crosses the atomic medium, we can also use quantum scattering theory. The 
Hamiltonian of the problem is 

H = Ho + V , Ho = Ha + Hp, (25) 

and we consider the initial state where all atoms are in their ground state and where a 
single photon of wave vector k = kuz and polarisation e = lia; is incident on the atomic 
medium 

\^,) = \g)®\k,e), (26) 
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with \Q) = \1 : g, 2 : g, ...,N:g). The state is an eigenstate of Hq with 
energy hu. The interaction of the photon with the atomic medium, described by the 
couphng V, can be viewed as a colhsion process during which an arbitrary number of 
elementary scattering events can take place. Each event starts from a state \Q) ® \q, s) 
and corresponds to: 

(i) The absorption of the photon in mode q, s by atom j, which jumps from its ground 
state \j : g) to one of its excited states \j : Ca)- The state of the system is then 

\£j,a) = \l-g, ...,j:ea,...,N:g)(S)\vac), (27) 

where |vac) stands for the vacuum state of the electromagnetic field. The subspace 
spanned by the states \Sj^a) has dimension 3A^. 

(ii) The emission of a photon in the mode {q' , s') by atom j, which falls back into its 
ground state. 

Finally a photon emerges from the atomic sample, and we want to determine the 
probability amplitude to find this photon in the same mode \k, e) as the initial one. 
The T matrix defined as 

where 0+ is a small positive number that tends to zero at the end of the calculation, 
provides a convenient tool to calculate this probability amplitude. Generally 

T,f = {^j\TiE,)\^,) (29) 

gives the probability amplitude to find the system in the final state |\E'/) after the 
scattering process. The states |\Efj) and are eigenstates of the unperturbed 

Hamiltonian Hq, with energy Ei. Here we are interested in the element Tu of the T 
matrix, corresponding to the choice |\E'/) = {"^i)- Using the definition (28) we find 

T---^j:^-''''-'HerAj—^je,^.). (30) 

jj' 

We now have to calculate the (3iV) x (3A^) matrix elements of the operator 1 / (z—H), 
with z = hw + iO+, entering into ( [30| ). We introduce the two orthogonal projectors P 
and Q, where P projects on the subspace with zero photon, and Q projects on the 
orthogonal subspace. We thus have 

P\£j,^) = \£j,^) P\g)®\k,e) = 0, (31) 

Q\£j,a) = Q\g)^\k,e) = \g) e). (32) 

We define the displacement operator 

and use the general result [20] 
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where the effective Hamiltonian ifefr is 

H,s = PiHo + R{z))P. (35) 

For the following calculations, it is convenient to introduce the dimensionless matrix 
[M] proportional to the denominator of the right hand side of (34): 

2 



[M]( 



Hcs\£j,a)- 



(36) 



J(i',«'),0» - 

It is straightforward to check[|] that for z ^ hu this matrix coincides with the symmetric 
matrix appearing in ( [l6) ). Indeed the matrix elements of R{z) are 



{£j',a'\Riz)\£j,a) 



q,s 



Cq S^Sa' 



(37) 



which can be calculated explicitly. For j = j', the real part of this expression is the 
Lamb shift that we reincorporate in the definition of Uq, and its imaginary part reads: 



KT 

{£j,a'\R{z)\£j^a) = -i^S, 



(3J 



For j 7^ j' , the sum over {q, s) appearing in (37) is the propagator of a photon from an 
atom in Vj in internal state \ea), to another atom in Vj' in internal state \ea')- This is 
nothing but (up to a multiplicative coefficient) the expression that we already introduced 
for the field radiated in r^/ by a dipole located in r^: 



{£j>^a'\R{z)\£j,a) 



-9c 



(39) 



where the tensor is defined in Eqs. (13 14). 

Suppose now that the atoms are uniformly distributed over the transverse area L^Ly 
of the quantisation volume. We set n^^"^^ = N/^L^Ly) and we rewrite the expression 
(30) of the desired matrix element Ta as 
TiiL, 1 



he 



2N 



(col) 



k{z 



1 



-don 



(COI)-Q 



(40) 



where the coefficient 11 has been defined in (22). The result (40) combined with (21) 
leads to 



r = 1 



he 



(41) 



which constitutes the 'optical theorem' for our slab geometry, since it relates the 
attenuation of the probe beam T to the forward scattering amplitude Ta. 

The emergence of resonant van der Waals interactions is straightforward in this 
approach. Let us consider for simplicity the case where only N = 2 atoms are present. 
The effective Hamiltonian ifgfr is a 6 x 6 matrix that can be easily diagonalized and 
its eigenvectors, with one atom in |e) and one in \g), form in this particular case an 
orthogonal basis, although i^eff is non-Hermitian [521133]. For a short distance r between 



X As for the derivation leading from Eq. ( 10 ) to Eq. ( 12 ), one must take into account the non-resonant 
terms that are usually dropped in the RWA, in order to ensure the proper convergence of the sum (37) 
and obtain the tensor gaa'- 
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the atoms {kr <^ 1), the leading term in hi{u) and h2{u) is and the energies (real 
parts of the eigenvalues) of the six eigenstates vary as ~ ±hr/{kr)^ (resonant dipole- 
dipole interaction). The imaginary parts of the eigenvalues, which give the inverse of 
the radiative lifetime of the states, tend either to F or when r — )• 0, which correspond 
to the superradiant and subradiant states for a pair of atoms, respectively [5^ . 

For N > 2 the eigenvectors of the non-Hermitian Euclidean matrix Hgg are in 
general non orthogonal, which complicates the use of standard techniques of spectral 
theory in this context [29l l30] . More precisely, one could think of solving the linear 



system (16), or equivalently calculating Ta in Eq. (30), by using the expansion of the 



column vector \Y) defined in Eq. (17) on the left (|aj)) and right ((/3j|) eigenvectors 
of -f^eff- Then one could inject this expansion in the general expression of the matrix 
element Ta, to express it as a sum of the contributions of the various eigenvalues of 
Hes- However the physical discussion based on this approach is made difficult by the 
fact that since H^s is non-Hermitian, the and the bases do not coincide. 

Hence the weight {/3j\Y){Y\aj) of a given eigenvalue in the sum providing the value of 
Ta is not a positive number, and this complicates the interpretation of the result. 



3.4- Beyond the sparse sample case: 3D vs. 2D 

For a sparse sample, we already calculated the optical density at first order in density 
(Eq. (24)) and the result is identical for a strictly 2D gas and a thick one. The 
approach based on quantum scattering theory is well suited to go beyond this first 
order approximation and look for differences between the 2D and 3D cases. The basis 



of the calculation is the series expansion of Eq. ( 34 ) , which gives 

P 



P- 



1 



z-H 



:P 



Hn 



n=l 



P 



Z-Hn 



PR{z)P- 



1 



Hn 



(42) 



Consider the case of a resonant probe 6 = for simplicity. The result T ~ 1 — (To^'-^'V^ 



obtained for a sparse sample in Eq. (24) corresponds to the first term [P/{z — Hq)] of 
this expansion. Here we investigate the next order term and explain why one can still 
recover the Beer-Lambert law for a thick (3D) gas, but not for a 2D sample. 



Double scattering diagrams for a thick sample (ki ^ Ij. We start our study by adding 
the first term [n = 1) in the expansion (42) to the zero-th order term already taken 



into account in Eq. (24). This amounts to take into account the diagrams where the 



incident photon is scattered on a single atom, and those where the photon 'bounces' 
on two atoms before leaving the atomic sample. Injecting the first two terms of the 



expansion (42) into (40), we obtain 
TuL-, 1 



he 



-aon 



(col) 



-t + 



N 



EE' 



,ik{zj-z-,) 



9x 



'33 



(43) 



We now have to average this result on the positions of the atoms j and j' . There 
are N{N — 1) A^^ couples (j, j')- Assuming that the gas is dilute so that the average 
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distance between two atoms (in particular \zj — Zj/\) is much larger than k ^, the leading 



term in g^^ is the e*"/M contribution of hi{u) in Eqs. (13)-(14). We thus arrive at 
TiiLz 



he 



-aon 



(col) 



3A^ 



ik{z-z')'_ 



^ik\r~r'\ 



(44) 



where the average is taken over the positions r and r' of two atoms. We first calculate 



the average over the xy coordinates and we get (cf. Eq. (20)) 
TjiL^_ 1 



he 



-aoTi 



(col) 



-aon 



(col) / ik(z—z') ik\z—z'\ 



(45) 



For a thick gas {ki ^ 1) the bracket in this expression has an average value of ~ 1/2. 
Indeed the function to be averaged is equal to 1 if 2; < z', which occurs in half of the 
cases, and it oscillates and averages to zero in the other half of the cases, where z > z'. 
We thus obtain the approximate value of the transmission coefficient: 

.TiiL-,. ^ 1 /-^^n 1 / r„^n\2 



fc£>l: r=l-i'- 



he 



Z o 



(46) 



where we recognize the first three terms of the power series expansion of T 
exp(— cro^'''^°^V2); corresponding to the optical density V = aQn^'^°^\ 



Double scattering diagrams for a 2D gas (i = 0). When all atoms are sitting in the 
same plane, the evaluation of the second order term (and the subsequent ones) in the 
expansion of Ta in powers of the density is modified with respect to the 3D case. The 



calculation starts as above and the second term in the bracket of Eq. (43) can now be 
written 



'J J 



n 



(2D) 



(47) 



If we keep only the terms varying as e^^/u in hi and /i2 (Eq. (14)), we can calculate 
analytically the integral in (47) and find the same result as in 3D, i.e., iaon^'^^^ A. If this 



was the only contribution to (47), it would lead to the Beer-Lambert law also in 2D, at 



least at second order in density. However one can check that a significant contribution to 
the integral in (47) comes from the region u = kr < 1. In this region, it is not legitimate 
to keep only the term in e^^^ /kr in hi, h2, since the terms in e^^^/{krY, corresponding 
to the short range resonant van der Waals interaction, are actually dominant. Therefore 



the expansion of the transmission coefficient T in powers of the density differs from (46 ) 



and one cannot recover the Beer-Lambert law at second order in density. Calculating 
analytically corrections to this law could be done following the procedure of [12|. Here 
we will use a numerical method to determine the deviation with respect to the Beer- 



Lambert law (see section 4.2). 



Remark. For a 3D gas there are also corrections to the second term in Eq. (45) due 



the 1/r^ contributions to hi and /i2. However these corrections have a different scaling 
with the density and can be made negligible. More precisely their order of magnitude 
is ~ n^^^^k~^, to be compared with the value ~ n^^°^^k~'^ of the second term in Eq. (45). 
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Therefore one can have simuhaneously n^^^^k ^ ^ 1 and n^^°^^k ^ > 1, if the thickness 
£ of the gas along z is 

4. Absorption of light by a slab of atoms 

In order to study quantitatively the optical response of a quasi-2D gas, we have 



performed a Monte Carlo calculation of the transmission factor T given in Eq. (21), 
and of the related optical density V = In |T|~^. We start our calculation by randomly 
drawing the positions of the N atoms, we then solve numerically the 3A^ x 3A^ linear 



system (16), and finally inject the result for the N dipoles in the expression of T. 

The atoms are uniformly distributed in a cylinder of axis z, with a radius R and a 
thickness £. The largest spatial densities considered in this work correspond to a mean 
inter-particle distance ~ k~^. Around each atom we choose a small excluded volume 
with a linear size a = 0.01 k~^. We varied a by a factor 10 around this value and checked 
that our results were essentially unchanged. Apart from this excluded volume we do 
not include any correlation between the positions of the atoms. This choice is justified 
physically by the fact that, in the case of large phase space densities which motivates our 
study, the density fluctuations in a 2D Bose gas are strongly reduced and the two-body 
correlation function g2{r, r') is such that (72(^7 t) ~ 1 [35j. 

In this section we first determine the value of N that is needed to reach the 
'thermodynamic limit' for our problem: for a given thickness £, V should not be an 
independent function of the number of atoms N and the disk radius but should 
depend only of the ratio N/irR^ = n^'^'^^\ We will see that this imposes to use relatively 
large number of atoms, typically N > 1000, for the largest spatial densities considered 
here. All subsequent calculations are performed with N = 2048. We then study the 
dependence of V with the various parameters of the problem: the column density n^^"^\ 
the thickness of the gas i, and the detuning A. In particular we show that for a given 
^(coi) recover the 3D result nl\ when the thickness i is chosen sufficiently large. 



4.I. Reaching the 'thermodynamic limit' 

We start our study by testing the minimal atom number that is necessary to obtain a 
faithful estimate of the optical density. We choose a given value of = N/nR^ and 
we investigate how V depends on N either for a strictly 2D gas {£ = 0) or for a gas 
extending significantly along the third direction {£ = 20 k~^). We consider a resonant 
probe for this study (A = 0). We vary N by multiplicative steps of 2, from A = 8 up 
to A = 2048 and we determine how large A must be so that V is a function of n*^™'-* 
only. 

The results are shown in Fig.|2^ and Fig.[2]3, where we plot P as a function of 
A. We perform this study for four values of the density n''™'\ corresponding to 
(Jq7^(co1) _ 2 and 4. Let us consider first the smallest value aQn^^°^^ = 0.5. For 
each value of A we perform a number of draws that is sufficient to bring the standard 
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Figure 2. Variation of the optical density T) — \n |T| 

function of the number of atoms N, for £ = (a) and £ — 20k~^ (b), and for 4 values 
of the density: aon'-^^^ ~ 0.5 (black), 1 (red), 2 (green) and 4 (blue). The bars indicate 
the standard deviations. The dotted lines give the value obtained for our largest value 
of N {N — 2048). The results have been obtained at resonance (A = 0). 



error below 2 x 10~^ and we find that the calculated optical density is independent of 
(within standard error) already for > 100, for both values of i. Consider now our 
largest value a^n^^"^^ = 4; for a strictly 2D gas {i = 0), V reaches an approximately 
constant value independent of N for > 1000. For aon^^°^^ = 4 and a relatively 
thick gas {i = 20 k~^, blue squares in Fig.|2]3), reaching the thermodynamic limit is 
more problematic since there is still a clear difference between the results obtained with 
1024 and 2048 atoms. This situation thus corresponds to the limit of validity of our 
numerical results. In the remaining part of the paper we will show only results obtained 
with = 2048 atoms for column densities not exceeding aon'^™'^ = 4. The number of 
independent draws of the atomic positions (at least 8) is chosen such that the standard 
error for each data point is below 2%. 

4-2. Measured optical density vs. Beer-Lambert prediction 

We now investigate the variation of the optical density V = ln|T|^^ as function of 
the column density of the sample 

^(coi)^ or equivalently of the Beer-Lambert prediction 
T^BL = n^^"^^cr. We suppose in this section that the probe beam is resonant (A = 0), 
and we address the cases of a strictly 2D gas {£ = 0) and a thick slab [i = 20 k~^). 

Consider first the case of a strictly 2D case, i = 0, leading to the results shown in 
Fig.|3|i. We see that V differs significantly (~ 25%) from V^l already for Dbl around 
1. A quadratic fit to the calculated variation of V for aon^"^^^ < 1 (continuous red line) 
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Figure 3. Variations of the optical density T) as function of the Beer-Lambert 
prediction 2?bl for ^ = (a) and £ = 20 k^^ (b). The black dotted line is the straight 
line of slope 1. In (a) the continuous red line is a quadratic fit 2? = Vbl (1 — M^^bl) 
with /i = 0.22 to the data points with I?bl < 1- The calculations are done for 
N — 2048, A = and the bars indicate standard deviations. 



gives 

Pf^^I?BL(l-0.22pBL). (48) 

The discrepancy between V and Pbl increases when the density increases: for "Dbl = 4, 
the calculated V is only ^ 1.4. For such a large density the average distance between 
nearest neighbours is ~ and the energy shifts due to the dipole-dipole interactions 
are comparable to or larger than the linewidth F. The atomic medium is then much 
less opaque to a resonant probe beam than in the absence of dipole-dipole coupling. 

Consider now the case of a thick sample, i = 20k~^ (Fig.jsjD). The calculated 
optical density is then very close to the Beer-Lambert prediction over the whole range 
that we studied. This means that in our chosen range of optical densities, the mean-field 
approximation leading to Vbl is satisfactory as soon as the sample thickness exceeds a 
few optical wavelengths A = 271 /k. 

It is interesting to characterize how the optical density evolves from the value for a 
strictly 2D gas to the expected value from the Beer-Lambert law Pbl when the thickness 
of the gas increases. We show in Fig. |4] the variation of V as function of i for three 
values of the column density corresponding to Pbl = 1,2 and 4. An exponential fit 
V = a + (3 exp{—i/ic) to these data for 2 k~^ < i < 20 k~^ gives a good account of the 
observed variation over this range, and it provides the characteristic thickness ic needed 
to recover the Beer-Lambert law. We find that ic ~ 3.0 k~^ for Dbl = 1, ~ 3.5 A;~^ 
for Vbl = 2, and 4 ~ 4.4 k'^ for Vbl = 4. 
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Figure 4. Variation of T) with the thickness i of the gas for various column densities 
corresponding to I?bl = 1 (red), 2 (green), 4 (blue). The continuous lines are 
exponential fits to the data. The dotted lines give the Beer-Lambert result. The 
calculations are done for N = 2048, A = and the bars indicate standard deviations. 



Remark. For the largest value of the column density considered here 



in 



(col) 



4) 



we find that V increases slightly above the value Pbl when i is chosen larger than 
20 (upper value considered in Fig. ffl). We believe that this is a consequence of the 



edge terms that we neglected when approximating Eq. (19) by Eq. (21). These terms 
become significant for Dbl = 4 because for our atom number = 2048, the sample 
radius R ~ 55 is then not very large compared to its thickness for i > 20 k~^. In 
order to check this assumption, we also calculated numerically the result of Eq. (19) 
(instead of Eq. (21)) for practical values of the parameters (position and radius) of the 
lens represented in Fig. [T^. The results give again I) I^bl, but now with V remaining 
below Vbl- Since our emphasis in this paper is rather put on the 2D case, we will not 
explore this aspect further here. 



4-3. Absorption line shape 

Resonant van der Waals interactions manifest themselves not only in the reduction of the 
optical density at resonance but also in the overall line shape of the absorption profile. 
To investigate this problem we have studied the variations of V with the detuning of the 
probe laser. We show in Fig. [s^ the results for a strictly 2D gas {i = 0) for n'^^^Vo = 1,2 
and 4. Several features show up in this series of plots. First we note a blue shift of the 
resonance, which increases with n^"^^^ and reaches A ^ r/4 for aon^'^^^ = 4. We also 
note a slight broadening of the central part of the absorption line, since the full-width 
at half maximum, which is equal to F for an isolated atom, is ~ l.SF for 
Finally we note the emergence of large, non-symmetric wings in the absorption profile. 
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Figure 5. (a) Variation of D with the reduced detuning 6 = 2A/T of the probe laser 
in the case of a 2D gas {£ = 0), for three values of the column density aon^'^°^'> = 1 
(red), 2 (green), and 4 (blue), (b) Blue full squares: same data as in (a), now plotted 
for 2?/2?BL as function of 6. Blue open squares; 2?/2?bl for a thick gas {£ = 20 k"^). 
Black stars: V/Vbl for a 2D gas {£ = 0) and a large exclusion region around each 
atom (a = fc~^). All data in (b) correspond to aon^'^^^ = 4. The calculations are done 
with N ~ 2048 atoms and the bars indicate standard deviations. 



This asymmetry is made more visible in Fig. |5)d, where we show with full blue squares 
the same data as in Fig. [s^i for n'^™'Vo = 4, but now plotting V/V^l as function of 
6. For a detuning 6 = ±15, the calculated optical density exceeds the Beer-Lambert 
prediction by a factor 4.1 (resp. 2.8) on the red (resp. blue) side. 

In order to get a better understanding of these various features, we give in Fig. [SJd 
two additional results. On the one hand we plot with empty blue squares the variations 
of P/Pbl for a thick gas {i = 20 k~^) with the same column density There 
are still some differences between V and Vbl in this already pointed out in 

[T7j . but they are much smaller than in the i = case. This indicates that the strong 
deviations with respect to the Beer-Lambert law that we observe in Fig. [5^ are specific 
2D features. On the other hand we plot with black stars the variations of P/Pbl for 
a 2D gas {£ = 0) in which we artificially increased the exclusion radius around each 
atom up to a = instead of a = 0.01 (blue full squares) for the other results in 
this paper. This procedure, which was suggested to us by Robin Kaiser, allows one to 
discriminate between effects due to isolated pairs of closely spaced atoms, and many- 
body features resulting from multiple scattering of photons among larger clusters of 
atoms. The comparison of the results obtained for a = 0.01 and a = k~^ suggests 
that the blue shift of the resonance line, which is present in both cases, is a many-body 
phenomenon, whereas the large amplitude wings with a blue-red asymmetry, which 
occurs only for a = 0.01 k~^, is rather an effect of close pairs. 
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This asymmetry in the wings of the absorption hne in a 2D gas can actually be 
understood in a semi-quantitative manner by a simple reasoning. We recall that for 
two atoms at a distance r ^ k~^, the levels involving one ground and one excited atom 
have an energy (real part of the eigenvalues of i^eff) that is displaced by ~ ±hT/{kr)^. 
A given detuning S can thus be associated to a distance r between the two members 
of a pair that will resonantly absorb the light. To be more specific let us consider a 
pair of atoms with kr <^ 1, and suppose for simplicity that it is aligned either along 
the polarization axis of the light (x) or perpendicularly to the axis (y). In both cases 
the excited state of the pair that is coupled to the laser is the symmetric combination 
(|1 : g;2 : Cx) + \l : e^; 2 : g))/\/2. If the pair is aligned along the x axis, this 
state has an energy hwo — 3hr/2{kr)^, hence it is resonant with red detuned light such 
that 6 = —3/{kr)^. If the pair axis is perpendicular to x, the state written above 
has an energy hujQ + 3hr/4{kr)^, hence it is resonant with blue detuned light such that 
6 = 3/2{kr)^. This clearly leads to an asymmetry between red and blue detuning; indeed 
the pair distance r needed for ensuring resonance for a given 5 > 0, ruue = (3/2|5|)^/^, is 
smaller than the value rred = (3/|5|)^^^ for the opposite value —5. Since the probability 
density for the pair distance is V{r) oc r in 2D for randomly drawn positions, we expect 
the absorption signal to be stronger for —6 than for +6. In a 3D geometry the variation 
of the probability density with r is even stronger {V{r) oc r^), but it is compensated by 
the fact that the probability of occurrence of pairs that are resonant with blue detuned 
light is dimensionally increased. For example in our simplified modelling where the pair 
axis is aligned with the references axes, a given pair will be resonant with blue detuned 
light in 2/3 of the cases (axis along y or z) and resonant with red detuned light only 
in 1/3 of the cases (axis along x). This explains why the asymmetry of the absorption 
profile is much reduced for a 3D gas in comparison to the 2D case. 

5. Summary 

We have presented in this paper a detailed analysis of the scattering of light by a 
disordered distribution of atoms in a quasi-two dimensional geometry. The particles 
were treated as fixed scatterers and their internal structure was modeled as a two-level 
system, with a J = ground state and a J = 1 excited state. In spite of these simplifying 
assumptions the general trend of our results is in good agreement with the experimental 
finding of [7j, where a variation of the measured optical density similar to that of Fig. |3] 
was measured. 

Several improvements in our modeling can be considered in order to reach a 
quantitative agreement with theory and experiment. The first one is to include the 
relatively complex atomic structure of the alkali-metal species used in practice, with 
a multiply degenerate ground state; this could be done following the fines of [211 122] . 
A second improvement consists in taking into account the atomic motion. This is in 
principle a formidable task, because it leads to a spectacular increase in the dimension 
of the relevant Hilbert space. This addition can however be performed in practice in 
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some limiting cases, for example if one assumes that the particles are tightly bound in a 
lattice [36|,|37]. When the atom-light interaction is used only to probe the spatial atomic 
distribution of the gas, neglecting the particle motion should not be a major problem. 
Indeed the duration of the light pulse is quite short (~ 10 microseconds only). Each 
atom scatters only a few photons in this time interval and its displacement is then smaller 
than the mean interatomic spacing for the spatial densities encountered in practice. The 
acceleration of the atoms under the effect of resonant van der Waals interaction should 
also have a minor effect under relevant experimental conditions. Finally another aspect 
that could be valuably studied is the interaction of the gas with an intense laser beam 
[55] . One could thus validate the intuitive idea that saturation phenomena reduce the 
effects of resonant van der Waals interactions [391 [8] , and are thus helpful to provide a 
faithful estimate of the atomic density from the light absorption signal. 
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